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abstract: 



If an elas-tfc bar vii'th srome ihitial- imperfection is subjected to 
increasing compressive axi'all bond,, the' lateral deflection, .measured 
an3Twiiere along its lenglh^. rncreases: monotonically; the ^ primary buck- 
ling mode becoming more and more: dominanl^ in the deflected shape. . A 
classical approximate technique,, due: -ter Ri . V. . Southwell^ .correlates 
axial load and lateral deflection for determination of : the primary 
buckling load of the. bar.. Hov;ever,. Southwell’s approximation may be 
inaccurate if the primary mode does, nodt: predominate in: the deflected 
shape. The present thesis proposes a: technique which employs axial 
loading^ but also transverse loading to insure dominance of the 
primary mode in the deflected shape. 

The technique allows experimental determination of the second 
buckling load, again using appropriate transverse loads to insure 
dominance of that mode in the deflected shape. 
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NOMEN GLATUEE: 



a 

n 


Coeffxcxents which: characterize deflection’ due to loads 


b 

n 


Coefficients which characterize initial -deflection 


e. 

1 


Error associated with the. ith data point* 


h 


ffeighf of column, cross— section (in y direction) 


i,j,n 


Inrteger subscripts 


k 


Lateral stiffness of column (lateral load to produce 
unit lateral deflection) 


V 


Deflection divided by corresponding load.(6/P) 


w 


Width of column cross-section 


X 


Axis parallel to applied axial load 


y 


Axis* parallel to applied, lateral load 


B 


Flexural rigidity El 


c 

n 


Compact notation for specified function of T and 


E 


Young ^s modulus 


I ,J ,K 
n n n 


Constants arising from evaluation of certain specified 
integrals 


L 


Effective length of column 


P 


Compressive axial load 


P 

n 


Buckling load for the nth mode 


Q 


Lateral load 


S 


Sum of the squares of e^ 


T 


Tensile axial load 


U 


Strain energy 


V 


Voltage 


w 


Work 


Y 

n 


Orthogonal functions 
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Measured^ defXectxan. 



1/k 
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I . INTEQD[JCTIO^r 



The problem of determining the elastic buckling* load of:^ a** column 
has been one of both theoretical interest: and practical importance.. 
Euler’s theory for initially straight columns has been, accepted and 
widely used for many years. R* V* Southwell L5-1 ‘treated the: problem' 
of initially imperfect columns from both theoretricaX and' experimental 
viewpoints. If, before loading, the column has some small distortion, 
usually a superposition of many modes, the applicatibn ol compressive 
axial load P magnifies the amplitude of the nth mode hy a factor 
P /(P -P), v/here P is the nth mode buckling load.. As: E increases,, 
magnification of the higher mode imperfections (n > 1) becomes* rela-- 
tively insignificant as compared to that of the first mode. If the 
column is prismatic and the deflection is measured at mid-span, even 
numbered modes make no contribution at all, and the dominance of the 
first mode is even more pronounced. Southwell postulated that to a 
close approximation <5 = 6^/(i-p/p^)^ where is the first mode imper- 
fection and <5 is the mid-span deflection. Alternatively P^v -6 = 6^ 
v/here v = 6/P. Accordingly, if a plot of v vs . 6 is made, a straight 
line with inverse slope equal to P^ results, provided P is large 
enough to assure first mode dominance. 

Successful use of Southv^ell’s method for determining buckling 
loads depends upon the unloaded column having an initial lateral 
deflection in which there is a significant component of the first 
buckling mode. If this is not the case, the desired linearity of 
(an appropriate) plot vj-lll not be observed and it will not be 
possible to make an accurate determination of the first buckling 
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load. This paper presents, an. extension ofr Southwell* s method:.. The 
proposed method does not reXy* on. accldeirtall imperfection. . Lateral. 
loads are used to insure dominance o£~ the nth mode in the deflected 
shape when determining the nth buckling load.. By using appropriate 
lateral loads, the buckling load for any mode may be determined,, 
although only the first two are considered:, here . 
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II. THEORETICAL ANALYSIS 




Figure 2 . 1 
Loading Diagram 



A. GENERAL CASE 

Figure 2.1 represents a bar, with either variable or constant 
cross-section 5 having some initial imperfection and subjected to 
axial loading P and lateral loading Q. All forces and axes are 
shown in their positive sense. The solution for the deflected shape 
may be found using a linear combination of orthogonal functions to 
represent the initial deflection and the deflection under load. The 
case represented here is for a single transverse load, but the method 
may be extended for any number of loads. Here and elsewhere it is 
assumed that all loads and deflections are confined to a single plane, 
namely the x-y plane. 

Let y be the initial ordinate, measured from the chord AB and let 
y^ be the deflection due to external loads. The total ordinate after 
the loads have been applied is y = y^ + y^. 

Assume y = b.Y. + b_Y_ + b.Y. + ... (2,1) 

■^o 1 1 2 2 3 3 

Yl = + a^Y^ + ... (2.2) 



11 



where Y = Y (x) (n = 1,2, 3:^^^..)' Is: a: srertr cdf ocrthogonal functions _ 
n n 

in the interval (0,L) whichL satisify the, dijEfenentlal. equation B-ia: 
subject to conditions B-lb (see Appendix B).. As" the loads are 
applied, the distance bet^^'een A and B decreases by an amount 



d(7i+yo) 



= jJ [Cy'f - (yp^J dx 



1 


rdy 1 


1 H 


dx. 



dx 



where y' = ^(a.+b.) Yl; Yi - ^a.Yl;: y^' — ^b.Y.l 

11 1^1 ir -^Q 3_1L 

and the prime denotes differentiation with, respect to x. 



(2:3) 



(y»)^ = EZ(a.tbJ(a.+b JYI: YV.;: (y’)^=^ ZIb .. h“.. Yi Yl 

1 



a/ 



= E (a + a.bj / (yp dx 



a.^ 



= E ( a + a.b.) J 
2 1 1 n 



(2.4) 



For definitions and properties of the functions appearing here and 
below, see Appendix B. 

Now consider the energy variations due to some small additional 

deflection da Y and assume da is positive in the same sense as Q. 
n n n ^ 

The work done by the force P is 
3(X- X ) 

dW^ = P — 5 a ' da = P(a +b )• J * da (2.5) 

P oa n n n n n 

n 

The work done by the force Q (located at x = ) is 

dW^ = QY (x^)- da (2.6) 

Q n 0 n 
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The increase in s trra fn energy is 



dU = 



au 

3a 



da = -r — 
n da 



Er(x) 
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da 



rr 



3a 



^ B(y")Zdx 



da 



where B = ElCx) = BCx) 



g 


1 r 


'• r 


dU ^ , 


~ LZ a., a.. / 


B YV. Y' . dx 


3a 

-n 


2 L X jJ 


I •' 



n 



3a 



n 



1 2 

1 E a.. K., 

2 XX. 



da 



a. 



= a K da 
n n a 



(2.7) 



Equating the increase in strain, energy and the work done by the 



external forces , an equation for the coefficients a^ may be found 



P(a +b )J • da + QY (x^) • da 
n n n n n 0 n 



a = 
n 



QY (x.) + Pb J 
n 0 n n 



K 



PJ 



a K da 
n n n 



or since K /J = P 
n n n 



n 



Y (x„) + |- b 
K n b P n 
n n 

n „ 



( 2 . 8 ) 



y = Yi + y. 



z 



K Y 
n n 

1 



(x^) + b 
b n 

_ ^ 

P 



Y (x) 
n 



(2.9) 



n 

These terms may be ordered so that ^2 ^3 * * " 

It is obvious, for values of P approaching F^, that if Q is sufficiently 
large, the first term of the series is strongly dominant. Southwell, 
rather than applying a lateral load Q, depended on some accidental 
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imperfection such that h^ f However,, if th& nature of- the . initial, 
imperfection is such thal >» b^ for same: n*. » 1,. the firsti term~does : 
not dominate unless P is very* nearly- equal, toe F^.. Consequently, . the 
linear relation which Southwell’s analysis depends on may not* appear 
in plots of experimental data^ 

Equation (2.9) also demonstraties something: which does not* seem* to 
have been treated in previous literature,, n the: value ofithe^ applied 
lateral load Q is sufficientTy the first: terra will dominate even .if 
is not the largest and even, if P does nol approach P^. Furthermore, 
if the lateral load is applied near the calumn midlength, then 
Yi(Xo)/Ki >> Y^(Xq)/K 2 - Thus,, in effect,, the first mode imperfection 
may be increased by the use of lateral loading, , insuring dominance of:' 
the first mode in the deflected shape,, regardless of the value of. 

P(P ^ increasing the effective firsl mode imperfection to 

bf* = b^ + (QY^(Xq)/K^) , the accuracy and reliability of Southwell’s 
method can be greatly enhanced. 

Another experimental procedure, which involves observing deflection 
as a function of lateral load Q for a few particular values of P, will 
be described subsequently. 

For the special case of columns which are symmetric with respect to 
a normal plane through the center,, if the lateral load is applied at 
mid-span, non-zero values of the coefficients a^ occur only for 
n = 1,3,5, ... These coefficients attenuate rapidly (l>>a^/aj^>>a^/a^ . 
and it is evident that the first term of equation (2.9) dominates the 
others even if 1-P/P^ is not small. 

This suggests that the value of might be accurately approximated 
even if tensile loadings (i.e., negative values of P) are used. First 
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term dominance is diminished slightly, birt not sign lli'canlrXy ^ . foi: 
tensile loading. In the experimental program to ha dasxrrihed,. teixsile: 
loading was used partly for convenience, hut largely ta demonstrate 
what appears not to have been done before, namely to evaluate a 
compressive critical load without ever applying compression. 

Assuring that the applied axial compressive force F daes not closely 
approach the first buckling load permits the experimerrtalL determi-- 
nation of the second critical buckling load!! lo do this' it is necessary 
to assure dominance of the second term of equation (2. .9).. This can be 
done, again using lateral loads to ’^excite” mainly the second mode.. 
Generally, it is difficult to develop an experimental technique, which 
will permit accurate establishment of the true first mode shape: 
and also one which is orthogonal to it. However,, if the column is 
symmetric with respect to a central normal plane, clearly the mode 
shapes are alternately symmetric and antisymmetric. Thus, antisym- 
metrical loading can ^’excite" only the second and higher antisym- 
metrical modes. If the lateral loading consists of two oppositely 
directed forces near the quarter points, fourth mode excitation will 
be small and the second should dominate. 

B. PRISMATIC BAR 

Although there is no necessity for determining the buckling load 

of an elastic, prismatic bar by experimental means since the theory 

provides complete information for that case, it is still reasonable 

to select such a bar for experimental confirmation of the general 

analysis. The particular form of the preceding results for the case 

of a prismatic bar may *be easily found. The orthogonal functions are 

Y = sin mrx/L. The constants K and J are then 
n n n 
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(2:iQ) 



K = 

EL 






2L 



T _ ■n^TT'^ 

■ ^rT 



Equation (2.9) is then 

2Ql3 



nrrx. 



y = Yl + 7o = E 



Ern‘+7r‘^ 



sin 



-r^ + b 
L n 



Slit.- 



mrx 



I - 






( 2 . 11 ) 



(2:12) 



For first mode determination-y the load- Q coul'd' he. applied anywhere , in 
(0,L) and the deflection could be measured at the same or any other 
position in C0,L}. However, it is convenient to. do both at the. center, . 
X = L/2. If this is done, then, x = = L/2. and equation (2.12) becomes 



y = 



2QL 



EiTt 






+ • • • 



1 - 



1 -. 



9P, 



(2.13) 



If instead of a compressive load P, there is a tensile load T, clearly 
T = -P and equation (2.13) becomes 

1 \ 



y = 



_ Q 



+ • • • 



1 + ^ 81(1 
1 



9P. 



) 



t y- 



where the quantity k represents lateral stiffness 



k = ZI-n^/2L^ 



and 



y = 



3 + 



il + 



1 + 



9P. 



1 + 



25P, 



(2.14) 
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With the application of Q at mid-span, the bar will deflect essentially 
in a first mode configuration, higher mode contributions being small. 



C. ESTABLISHICNT OF REFERENCE LEVEL 



The last term in equation (2.14) presents a problem in that its 
exact value is unknown and cannot be determined in the laboratory. The 
problem may be circumvented however. Referring to Fig. (2.2), three 
load cases are shown: 

(a) T = Q = 0; 
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(c) 



Figure (2.2) 

Experimental Loading Cases 



(b) Q = 0, T = T; and (c) Q = Q, T = T. With no loads applied there 
is some initial imperfection b. When load T is applied (Q = 0), the 
amplitude of the imperfection is reduced to an amount y* which is 
exactly equal to the last term in equation (2.14). With the load T 
applied and maintained, a reference level may be established for all 
subsequent deflections due to lateral loading. Now as the lateral 
load Q is applied (holding T constant), the deflection measured in 
the laboratory is simply 






T T 

1 + — 81(1 + y ) 
1 1 



625(1 + 



25P, 



(2.15) 
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D. NON-PRISMATIC BARS 



For non-prisTnatic harSy equa^xon (2.31) dasrcrlbes the deflected 
shape in terms of the applied loads,, the iiriflall imperfection^ .and : 
some possibly unknown ortrh.agonal functrions Y^.. The forms of . 

Y^, ... are clearly not harmonic in. x,. but: it should still be possible 
to excite mainly the first chanactenistic mode Y^ using a lateral! 
load Q applied near the column mxdlength.. The_ deflected shape vwill_ 
not be a half-sine curve as for prismatic: bars:,, but it still!. v/illlnot' 
have any intermediate nodes.: 

If the reference level for laboratory* measurements is established 
as before and the measurements are taken at: the: point of load lapplica- - 
tion X = x ^5 the measured deflection, is then 



6 = 



QY^(Xq) I 









. 1 

Y^(Xo>K3 1 + ^ 



(2.16) 



E. ANALYSIS OF EXPERIMENTAL DATA 

Continuing with the general case of a non- prismatic bar, once 
experimental data has been obtained, it may be analyzed, using equation 

(2.16), by several methods. Two specific methods are outlined here. 

Considering the relative magnitudes of the series terms in equation 

(2.16), the first is clearly the dominant term for any given value of T. 
If only the first term is retained as an approximation, 6 is a linear 
function of Q for any given T since K^, Y^(x^) and P^ are all constant. 

A plot of 6 vs. Q, for a given T yields a straight line v;ith an inverse 
slope equal to 
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Thus kp is a linear fuiicrtiaii af ilrdicralriTr^ t±iat: a plot* off vs. . TT 

will have a slope equal to k/P^^ and an. ihtrercBpIr: equal to k.. By 

simply dividing the intercept by the slope a value of may be found. 

A second method, which may offer greatrer accuracy, is a least 

squares approach using P^ as a parameter.. Consider again equation 

(2.16). For beams which are ’^nearly prismalic” the functions Y^(x) 

will be "nearly^^ sine functions so that an. estimate oP' the relative- 

values of the first several coefficients Yf ^(x^)K^ /Y^ )K can: be 

n Oil On: 

made by assuming the functions are indeed, sines . Similarly the 

ratios P^/P^, estimated to be 4,. 9,. 25, 

This provides a way of estimating: the contributions of terms: beyond _the 
first. If these added terms can be at least estimated, they can' be 
retained and used advantageously;, if not,, only the first term is' 
retained. 

For any given value of the parameter P^, equation (2.16) maybe 
written 



k 






1 , 

\ 1 + 



— C. = Q. p C. 
k 1 11 



1 ^ 

1 , ^ " vJ(x„)K3 
^2 



T. 

3 



(2.18) 



where p = 1/k. In the series = C(T^), only those terms are retained 

which can be estimated with some confidence. Let e. be the error 

1 



associated with the ith data point. 



e. = Q. p C. 
1 1^1 



(2.19) 



Let S equal the sum of the squares of e^ for n data points. 



S = E ef 
i=l ^ 



( 2 . 20 ) 



19 



To find the most probable (in the sense of least:, squanes) value; of rp 



and 



from the experimental data^ proceed ats’ follows: 



1 3S 

2 9p 



0 



n 

2 Q. C. 
i=l 1 1 



e. 

1 



( 2 . 21 ) 



n 

Z 

i=l 



(Q^ p 



n 



Z 

i=l 



6. Q. C. 
11 1 



= 0 



n 

1 5. Q. C. 

i=l 1 1 1 



z (Q. c.y^ 
i=i ^ 1 



( 2122 ) 



(2123) 



Now find p and S for different values of tha parameter P^. The. quantity 
S is a minimum when the most probable value of is‘ used also" yield- - 
ing a most probable value of p. 

For the case of a prismatic bar,, equation (2.15) describes the 
deflected shape. The two methods just outlined are, of course, applicable 
to this simplified case. All of the terms in equation (2.15) are known 
and hence all may be used to improve the accuracy of the least squares 
analysis if desired. 

Both of the methods outlined above may be easily programmed for the 
digital computer. Appendix A gives program listings for the prismatic 
bar case. 



F. DETERMINATION OF THE SECOND BUCKLING LOAD 

Determination of the second mode buckling load using the theory 
outlined above is straightforward assuming that El(x) = El(L-x). If 
an upv^ard load -Q is applied at x = x^ and a doi'/m^xard load Q is applied 
at X = ^“^0’ resulting equation for deflection is similar to 

equation (2.12). (It is convenient to apply the lateral loads at 
X = L/4 and x = 3L/4.) Such lateral loading does not excite odd 
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numbered modes, onljr ewn numbeized: modes",, andi if: Q is sufficiently 
large, the term correspondijig to the second" mode can be made 'to 'dominate 

The first, third, tifthy terms are negligible, and the fourth, .sixth 

... are small compared to the second. Thus, .the sort of data' reduction 
discussed above will lead to an evalna'tion ol the second mode.. 

If El(x) ^ EI(L-x), that is if there, is:, no; symmetry with respect~to 
a normal plane at x = L/2y it may* be. quite' difficult to apply ailateral . 
loading system which does not excite the firslrmode. A suggested lap- - 
proach is to apply an upward load at x — 1/4- and an equal- dowmzard iload 
at a point ^tear” x = 31/4.. Unless the second point is selected: 
correctly, the mixture of first and second, modes thus excited, will, 
not provide a linear plot of data.. One would.'. have to try several- such 
points in order to obtain the desired sort of' plot which could then 
be interpreted so as to yield the second buckling load. 
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in. EXFERUmUTAL EUaCEDUKE: 



A. DESIGN 

During the design and fabrication, of tha tes±: apparatus , strict 
attention was paid to simplicity.. The major problems: which arose 
where how to best simulate pinned end conditions' irr. the laboratory, . 
and how to apply the required axial and laterai leads:.. Several.designs : 
were considered which could easily transmit an applied^ compressive . 
axial load, but they could not easily support an applied lateral load: 
without introducing unknown end moments.. Application of a compressive 
axial load using a universal testing machine: could- be accomplished, . 
but lateral loads could not be easily applied or measured. Dead weights 
seemed to be the simplest means of applying both axial and lateral 
loads. The final design of the test apparatus is shown schematically 
in Fig. (3.1a). 

B. FABRICATION 

Compressive axial loads could not be easily applied using dead 
weights so it was decided to use tensile loading for convenience. This 
also allowed determination of higher mode buckling loads. However, 
there is a penalty for using tensile axial loads. This will be discussed 
in Section IV- B. Pinned end conditions were simulated using steel 
cylinders on either end. The cylinders were polished and rested on 
a finished surface. A mixture of graphite, molybdenum sulfide, and 
machine oil was used as a lubricant to eliminate dry friction between 
these cylinders and the test bed. Flexible steel cable was used to 
transmit the axial load T. The cable v;as brazed to a steel ring for 
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attachment to the end cylinders. A slat: was cut: itr eacJr errdi cylinder,, 
perpendicular to the longitudinal axis and t±e rih^ v?an ihsBrted. . A 
pin which provided a ^'knife edge'*^ contad was insertred: paralTell to the- 
longitudinal axis of the cylinder, passing through the eye of:' Ihe 
ring, see Fig. (3.1b). The knife edge of the pin ran dire clly ■ along 
the longitudinal axis of the cylinder. The effedxve lengtrh of^ the* 
column was measured from centerline to centerline of the end cylinders, 
and was 19.00 inches. Since a near-zero moment exists near- these: 
centers, changes in El near the ends have negligible effect: on bending . 
and the effective column length may be safely considered as the dis.tance 
be-tween centers of the end cylinders . The column itself was 1/^"’ x 178^.' 
cold drawn 1018 steel, with a modulus of elasticity f — 29. .4 x ID.^ psil. 
It was brazed to the end cylinders. The lateral load was applied with 
a knife edge and dead weights as shown in the figure. Deflection 
measurements were taken with a simple wire probe attached to a verti- 
cally mobile micrometer. The knife edge assembly used for application 
of lateral loads had a hole drilled in the flat upper surface which was 
filled with mercury. The probe was lowered until contact with the 
mercury was made, completing a simple neon light bulb circuit, see 
Fig. (3.2). For second mode buckling tests, two lateral loads were 
applied, one at each quarter point. The downward load was applied as 
before. The upward load was also applied using a knife edge and dead 
weights. A pulley was used to convert dead weight into upward force. 
Since an upward force was applied at one end of the beam, the cylinder 
on that end had to be restrained against upward motion. This v/as done 
with a simple ”hood^\ see Fig. (3.3). 
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Figure 3.1b 

TOP VIEW of’ END CYLINDER 







APPARATUS FOR 



Figure 3.2 

DETERMNATION OF FIRST BUUKLING LOAD 
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Figure 3*3 

APPARATUS FOR DETERfUNATION OF SEGOI® BUGKLira LOAD 
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C. TESTS 



The reference level for testing was established- before each test 
by first applying the axial load T, and then taking a reading from 
the micrometer with the knife edge assembly in place. The weight 
of the knife edge caused some small deflection,, but this may be 
treated in the same manner as the initial imperfections were. Once 
the reference level was established, the lateral load was increased 
from 0.0 to 4.0 pounds in half-pound increments with a final 1/4 
pound increment to give a maximum load of 4..25 pounds. No greater 
lateral load was applied so as to assure that plastic bending did 
not take place. This was done for 15 different values of T ranging 
from 0.0 to 42.0 pounds, in three-pound increments. For each such 
loading, the values of loads and corresponding micrometer readings 
were recorded. 

The same method of establishing a reference was used for second 
mode buckling tests. An upward load (-Q) was applied at x = L/4 and 
a downward load (Q) was applied at x = 3L/4. The value of Q ranged 
from 0.0 to 6.5 pounds, in half-pound increments. The range of Q 
could have been extended another two pounds without causing plasticity, 
but the effects of friction became evident above 6.5 pounds, since 
increases in Q caused almost immeasurable changes in deflection. The 
values of T ranged from 0.0 to 41.0 pounds in five-pound increments, 
the odd pound arising from the weight pan used. A few anomalies in 
the data for second mode were noted and will be discussed subsequently. 
Several attempts to find the cause of the discrepancies were made, but 
proved fruitless . 



27 



Once the data was oh-ta.f ned.^. ilt: was: aaaiyzedz using both offthe: 



methods described previously .. 

As a means of substantiating earlier, test: results , another , type 
of experiment was done to determine the primary buckling load, using 
axial compression and interpreting the data: according to Southwellls 
method. The apparatus used for this test is^ shown: in Fig. (3.4). . A 
cantilever beam,, embedded in a. vertically * mobile: "wall" was used to: 
apply a compressive axiai load, to the column.. The cantilever was 
instrumented with four strain gages,, t/o. mounted parallel to the 
direction of principal strain,, and. two perpendicular to it to compen- 
sate for Poisson *’s effect.. The gages were, connected in a four-arm"' 
bridge circuit and hooked to a Baldwin Modell SR-4- Strain Indicator . -- 
Type N. The column was also instrumented . with strain gages. One gage 
was mounted on each side of the column at raid-span, parallel to the 
direction of principal strain. These two gages were designed specif- 
ically for use on mild steels and were self-compensating for small 
ambient temperature changes. These two gages were connected in a two- 
extemal-arm bridge hooked to an Ellis Bridge Amplifier and Meter 
(BAM-1). The output of the BAM was monitored on a digital voltmeter. 

The cantilever strain gages were calibrated using dead weight tip 
loads, providing a direct correlation between indicated strain and 
applied load. 

Before each test, the no-load reference for the cantilever bridge 
was noted. A small load was then applied to hold the column vertical, 
and the reference reading for the column bridge was noted. Preliminary 
tests indicated that the Southwell plot became linear for P .5P^. 
Therefore, the axial load, including the weight of the end cylinder, 
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was increased from 16. to 26. pounds in one-pound increments. At 
each load increment the strain readings from the cantilever bridge 
and the voltage readings from the column bridge were recorded, the 
strain reading being proportional to applied axial load, and the 
voltage change (AV = V - V^) being proportional to column deflection. 

AV/P versus AV was then plotted for each of four runs and the 
value of was determined by the reciprocal of the slope. A compatcr 
program was utilized to obtain the best fitted line, for each mm, 
again using a linear least squares fit, and the average value of P^^ 
for all runs was calculated (see Appendix A). 

A standard tensile 
used for the column to 
result was E = (29.4 + 



test was run on a specimen of the material 
determine Young’s modulus experimentally. The 
.1) X 10^ psi. 
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IV. RESULTS AtrD DISCUSS lOH 



A. RESULTS 

The figures tabulated in Table I represent: the experimentally 
determined values of the buckling load for both first: ancL sacond 
modes. As described earlier (secrtion 2-E) two methads were used: to: 
analyze the data: one being a Linear approximation,, hereafter' 

method I; the other being a least squares analysis retaiTxiTig' any 
desired number of terms in equation (2 ..15),, hereafter method. IT. . 

The average value of arrived at by use of Southwell* a method- is 
also tabulated. 



Southwell’s 

Method I Method II Method 

Number of 

terms retained 1 14 

1st mode 

buckling load (lbs.) 31.0 32.0 31.2 32.1* 

2nd mode 

buckling load (lbs.) 127.8 110.5 109.0 

Table I 

Summary of Results 

Typical plots of 6 vs. Q and k vs. T for method I are shown in Fig. 
(4.1) and (4.2) respectively. A typical plot for determination of 
by Southwell’s method is shown in Fig. (4.3). Compiled data for all 
tests is listed in TablesII - V. 

Average value of four tests. 
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DEFLECTION 5 (IN 




LATERAL LOAD Q (LBS.) 
Fig. (4.1) 
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Table II 
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2ND MODE BUCKLING DATA 
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DATA. FDR SaUTHWELE’^S- METEOE 
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AV/P 


CSR 


F 


V 


AV 


AV/P 


1226 


16.3 


- 56 


49 


3. .00 


1226 


16.3 


- 61 


53 


3.25 


1211 


17.3 


- 64 


57 


3.29 


1226 


16.3: 


- 66 


58 


3.35 


1197 


18.3 


- 72 


65 


3. .55 


1197 


18.3 


- 73 


65 


3.55 


1182 


19.3 


- 84 


77 


3..99 


1182 


19 . 3. 


- 81 


73 . 


3.78 


1168 


20.3 


- 94 


87 


4.29 


1168 


20.3 


- 92 


84 


4.14 


1153 


21.3 


-105 


98 


4- ..60 


1153 


21.3 


-106 


98 


4.60 


1139 


22.3 


-126 


119 


5.. 3 3 


1139. 


22.3 


-123 


115' 


5.16 


1124 


23.3 


-144 


137 


5..8S 


1124 


23 . 3- 


-140 


132: 


5.67 


1110 


24.3 


-163 


156 


6.42 


1110' 


24.3: 


-162 


154 


6.34 


1095 


25.3 


-186 


179 


7. .08 


1095 


25'. 3: 


-197 


191_ 


7755 


1081 


26.3 


-224 


217 


8.. 25 


1081 


26. 3- 


-227 


219 


8.33 




V = 
o 


-7tiiv 








V = 
o: 


-8mv 







CSR 


P 


V 


AV 


AV/P 


CSR 


F 


V 


AV 


AV/P 


1226 


16.3 


- 46 


45 


2..76 


1226 


16.3 


- 55 


50 


3.07 


1211 


17.3 


- 54 


55 


3.12 


1211 


17.3 


- 61 


56 


3.24 


1197 


18.3 


- 65 


64 


3.50 


1197 


18.3 


- 69 


64 


3.50 


1182 


19.3 


- 74 


73 


3.78 


1182 


19.3 


- 80 


75 


3.89 


1168 


20.3 


- 96 


95 


4.68 


1168 


20.3 


- 89 


84 


4.14 


1153 


21.3 


-106 


105 


4.93 


1153 


21.3 


-106 


101 


4.74 


1139 


22.3 


-116 


115 


5.16 


1138 


22.3 


-121 


116 


5.20 


1124 


23.3 


-143 


142 


6.09 


1124 


23.3 


-138 


133 


5.71 


1110 


24.3 


-159 


15 8 


6.50 


1110 


24.3 


-160 


155 


6.38 


1095 


25.3 


-180 


179 


7.08 


1095 


25.3 


-181 


176 


6.96 


1081 


26.3 


-209 


208 


7.91 


1081 


26.3 


-211 


206 


7.83 



V = “Imv V — -5mv 

o o 



Note; GSR = Cantilever strain reading (pin/in) 

P = Compressive axial load (lbs) 

V = Voltage reading from digital voltmeter (mv) 

An empirical equation based on earlier calibration tests was used to 
relate CSR and P 

P(lbs) = (CSR -CSR)xl.01b/14.76pin/intWt 
o c 

For all tests CSR = 1440pin/in. 

Wt^ = weight of tRe end cylinder =1.5 lbs. 

AV = V-V 

o 



Table IV 
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DATA FOR TEN’S ILE TEST TQ 
DETEKKINE YQDTTG’’S MODQLUS 



T 


SRIL 


SRDL 


T 


SIL 


SDL 


aspg: 


STRESS 


0 


1920 


1922 


0 


0 


0 


0 


o; 


100 


1700 


1700 


100 


220 


222 


no.. 6 


3226 


200 


1478 


1482 


200 


442 


440 


220.. 5. 


6452. 


300 


1255 


1262 


300 


665 


660 


33X..3: 


9677' 


400 


1038 


1043 


400 


882 


879 


440. .3T- 


12903; 


500 


819 


827 


500 


1101 


1095 


549.. o: 


16129; 


600 


602 


610 


600 


1318 


1312 


657.5 


19355; 


700 


382 


385 


700 


1538 


1537 


76 8. .8. 


22581 


800 


161 


163 


800 


1759 


1759 


87 9. .5 


25806 


850 


53 


53 


850 


1867 


1869 


934.. 0 


27419 


Note: 


Strain 


readings 


(uin/in) 


were taken 


with 


increasing- loacL 



(SRIL) and decreasing load (SRDL).. 



T = Tensile axial load (lbs) 

ASPG = Average strain per gage 

The cross-sectional area of the specimen was 0.031 in^. 

SIL = Actual strain, increasing load 

SDL = Actual strain, decreasing load 



Table V 



B. DISCUSSION 



In order to esIiaLIxsh. a. gcdhl:. of ramparasorr. between experimental _ 

and theoretical results,, the Eulerr load- may bev calculated, . based _on 

the following measured quantities : 

E - Young ^s Modulus ~ (29 .A t ..1) x 10^ prsi 
w = column width = Q..2M-E j- ..0!CQ1 in., 
h = column depth - Q..I25 t ..QOl in.. 

L = effective column. length. = 12.. O'O' t- .Dll in.. 



TT%r _ tt^E wh^ 



3-2. lbs.. 



(4.1) 



Due to the uncertainties in the measured: quantities , there is an'' 
uncertainty in the calculated, value crf'F^.. The magnitude of 'this 
uncertainty may be approximated using- a ^logarithmic" error analysis . 
Taking the logarithmic derivative of equation (4.1) 



1 — 1 
P-I 




dE 




dw 


+ 3 


dh 


+ 2 


dL 


p. 




E 




w 




h 




L 


1 



















dPi = P^(.0324) = 1.05 lbs. 

P^ = 32.44 ^ 1..05 lbs. (4.3) 

The experimental values of P^ as determined by methods I and II fall 
only slightly outside the range of uncertainty associated with the 
calculated Euler load. The value of P^ determined by Southwell *s method 
may be considered as an upper bound in view of Southwell *s theory [5] 
which indicates experimental values will be either exact or high. 

The Euler load for second mode buckling is four times that for 
first mode. Thus, we would expect 

P^ ^ 129.32 ^ 4.20 lbs. 
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The experimeatal value detrermined: by rneitrad: H is? well, within: the 
range of uncertainty hut: that: deteG^raihed. by ' method, II- seems con- - 
siderably in error^ Thin diacreqjancy ik attributable to inconsis-- 
tencies in measured deflection.. Referring;^ to' Tabie III, particularly 
the data for T = 11..Q and T ~ 1S..0 lbs..,, iir may be seen that increas- 
ing T resulted in increased deflection.. For every given value of:' 
lateral load Q, the deflections measured with T" ='16.0 lbs. are 
approximately ..Q03 indies greater than those: measured for T =“lli0' 
ibs. A second test with T — 16'..0' lbs., was run on another day as'a: 
check, and the values obtained earlier v;ere_- repeated. This discrepancy 
was completely unexpected in view ox the theory ,. so' several checks of. 
the experimental apparatus were made,, but noc unexpected sources of:' 
error came to light.. This inconsistency occurred for several (but 
not all) values of increased T. 

Despite these anomalies, the data can be treated using method I. 

For a given value of axial load (T), a plot of deflection (6) vs. 
lateral load (Q) is made but only the slope of this line is considered. 
For this reason a constant error in measured deflection has no effect 
on the outcome of the analysis.. Method 11,. on the other hand, relies 
on the absolute rather than the relative magnitudes of measured 
deflections, and gives erroneous values for as expected. 

The theory developed earlier, particularly equations (2.13) and 
(2.14), shows why it is advisable to retain several series terms when 
the axial loading is tensile. For any given tensile loading, the 
effect of additional series terms in equation (2.14) is greater than 
that for compressive loading (Eq. 2.13). Thus a single term approxi- 
mation may be more easily justified for compressive leads than for 
tensile. 
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CON'CLLTSrQNH AITO mCmi^mATTOliS^ 



Experimental determinaticm of the primary buckling load .of:' a: bar 
can be accomplished without reliance, on aome accidental, initial 
imperfection. Lateral loading is used to assure dominance of: the- 
primary mode in the deflecrted shape.. The axial load may be either* 
tensile or compressive* By ccjrrelating observed- deflectiDn and . 
applied lateral load for several distincrt values ol' axial load, .the 
buckling load may be easily and accurately determined.. 

If a tensile axial load is used,, experimental de termination:. of : 
the second buckling load is possible,, using appropriate lateral, loads 
to assure dominance of the second mode in the deflected shape. . 

The developed theory and experimental technique, are generally 
applicable to bars of any cross-section, although orthogonal load-- 
ing schemes for determination of higher mode buckling loads may be 
difficult experimentally for non-prismatic bars. The use of compress- 
ive axial loading is advantageous for non-prismatic bars since it 
minimizes the effects of possibly unknown higher modes in the deflected 
shape. 

It appears that techniques suggested may also be extended for use 
in determining plate buckling loads. 
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APPEirarX k - COMPUTES PROGRAM' LIETTRGS: 



The computer programs listed on the following pages v/ere used 
to analyze experimental data.. All programs, are written in PORTRAIT 
IV for use on an IBM 360 computer.. 

The programs for Method I and Southwell’s: Method are essentially 
the same. The program was taken, from McCalla: [33,- in part,. with, ae 
few additions and minor changes made as necessary.. 

The input required for each, program,, excluding data, is minimal _ 
and is explained in 
programs the symbol 



the body of each program.. (Note that in the 
P is used to denote axial- tension. ) 
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DO 22 1=1, NPl 
22 YC=VC+Y(I) 

FORM SUMS OF PRODUCTS Y(I)^X(D* 
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INVERT NORMAL MATRIX 



C IS NOW THE INVERSE OF C 

FORM MATRIX PRODUCT OF INVERSE AND RIGHT-SIDE MATRIX 
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»» t i i, 

viJF 




FIND THE ERROR ASSOCIATED WITH THE I TH DATA POINT 
SUMSQ IS THE SUM OF THE SQUARE OF THE ERRORS. 
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25 YX(K) =YX( K)+Y( I )*X( I )**K 
30 CONTINUE 

GENERATE NORMAL MATRIX C USING SUMS OF POWERS OF X(I) 
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APPENDIX B - evaluation: OF INTEGRALF DIVOLVING: OPTHOGONAL 



FOlfCTrONS Al^D THEIK DERIVATrVES; 



The equation and conditions governing- the, shape: oft the. strained 
centerline of a pin-ended column, subjecrted tot an axial compressive ■, 
load P are 

BY" + PY = Q where H = B(x)' = EI(x) (Brla). 



YCO) = Y(L) = 0 



(B-lb) 



Let (Y., P.) and (Y., P..) be two distinct solutions to this' equation 
1 1 D 1 



(Bf2) 



BY." 


+ P..Y.. 


= 0 


1 


1 X 




BY." 


+ P.Y. 


= 0 


I 


J I 





Multiplying both sides of equation B-2 by Y." 

L L ^ 



/ BY." Y." dx + 

■i ^ ^ 

0 0 

L L 



dx = 0 



/ P.Y.Y." 

J 1 1 3 
0 

L 

- P. / Y.Y." dx = / BY."Y." dx = -P. f Y."Y. dx 

° L ° 

I Y."Y. dx 

J ^ ^ 

-Jn 



0 

'L 



P. / Y.Y." dx = P. 

Ill I 



= P 



■'0 L 

i 

Y.Y. 

L 1. 



Y.'Y.' dx 
1 J 



= -P. / Y. 'Y,' dx 
3 1 ] 



Also the left side of equation B-4 is 

L L 

P. r Y.Y." dx = -P. r Y.'Y.' dx 

- V ^ ’ 

(P. - P.)/ Y. 'Y. ' dx = 0 
1 I J II 

0 



(B-3) 



(B-4) 
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But: P_ and F.. are distiircir- (ii P i') 

I- I 

L. 

•* dx = O' i- j 

Also, frorn equation B-4 
L L 

f Y\."Y .. dx - — / Y..!^Y... dx — 0 i 

J ^ 1 J ^ T 

0 0 

By equadxcn. E— 5 

L L 



r BY.."Y,"' dx = -R..r 

J ^ ^ 4 



dx = -P... Y...’^, dx = 0 i ^ i 

T X 



Q 0. 

Soimnarizing the preceedihg argtmien1:s 



(B-5) 



(B-6) 



(B-7) 



/■ T..Y',." dx - I..6 ;; E.. = Y.;Y." dx 

/ t I r X3 r 1 



1 



)^dx 



-dx 



/ Y.'Y.' dx = J.6. J. (Y.' 

4 / ' V " 

/ BY."Y."dx = K.6..; K. =/ (Y.") 

J X 3 ^ 



where 6 . . is the Kronecker delta 
13 



0 i j 

6 . . = ■ 

II 



Returning to equation B-3 and B-4 it may be seen that 



K. = P.J. 
1 11 
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